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The scaled particle theory is applied to a description of thermodynamic properties of one-dimensional hard rod 
fluid in disordered porous media. To this end, we extended the SPT2 approach, which had been developed pre- 
viously. Analytical expressions are obtained for the chemical potential and pressure of a hard-rod fluid in hard 
rod and overlapping hard rod matrices. A series of new approximations for SPT2 are proposed. It is shown that 
apart from two well known porosities such as geometrical porosity and specific probe particle porosity, a new 
type of porosity defined by the maximum value of packing fraction of fluid particles in porous medium should 
be taken into account. The grand canonical Monte-Carlo simulations are performed to verify the accuracy of 
the SPT2 approach in combination with the new approximations. It is observed that the theoretical description 
proposed in this study essentially improves the results up to the highest values of fluid densities. 
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1. Introduction 

It is a great pleasure for us to dedicate this paper to our good friend and colleague Orest Pizio - an 
excellent scientist in the study of fluids on complex surfaces. A typical example of such complex systems is 
a fluid confined in a disordered porous medium. Much theoretical effort has been devoted to the study of 
fluids in porous materials for the last two decades starting with a pioneering work by Madden and Glandt 
Q|. In this paper a simple model for a fluid adsorption in a disordered porous medium was proposed. 
Within this model, a porous medium is presented as a matrix composed of quenched configurations of 
model fluids. Using the replica Ornstein-Zernike (ROZ) integral equations \2], the statistical-mechanics 
approach of liquid state was extended to a description of different fluids confined in disordered porous 
matrices [ 3, 4]. An important contribution into the development of this approach and its application has 
been made by O. Pizio with co-authors I4ll3ll. Despite a comprehensive study of a fluid in disordered 
matrices, no analytical result has been obtained using the integral equation approach even for a simple 
model like a hard sphere fluid in a hard sphere matrix. The main complication in a description of such 
a model appears due to the absence of a direct interaction between particles from different replicated 
copies of a fluid. As a result, the description of this model is equivalent to the study of a hard sphere 
mixture with strongly non-additive diameters, for which it is very difficult to find a correct analytical 
result within the integral equation approach. 

The first rather accurate analytical results for a hard sphere fluid in hard sphere (HS) and overlap- 
ping hard sphere (OHS) matrices were obtained quite recently 1 14-17] by extending the scaled particle 
theory (SPT) flsWioll to a HS fluid confined in a disordered porous media. The SPT approach is based on 
a combination of the exact treatment of a point scaled particle in a HS fluid and the thermodynamical 
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consideration of a finite size scaled particle. Exact results for a point scaled particle in a HS fluid in disor- 
dered porous media were obtained in 1 14]. However, the proposed approach referred to as SPT1 has got a 
subtle inconsistency appearing when the size of matrix particles is essentially larger than the size of fluid 
particles. This inconsistency was eliminated and thus a new improved approach was developed referred 
to as SPT2 1 17]. The expressions obtained in SPT2 include two types of porosities. One of them is defined 
by the pure geometry of porous medium (geometrical porosity ^>o) and the second one is defined by the 
chemical potential of a fluid in the limit of infinite dilution (probe particle porosity (/>). Based on the SPT2 
approach, the SPT2b approximation was proposed that reproduces the computer simulation data with 
a good accuracy at small and intermediate fluid densities. However, the proposed approach encounters 
a serious problem at high densities of a fluid. The expressions obtained in SPT2b lead to a divergence 
when a packing fraction of a fluid reaches the value equal to the probe particle porosity. Consequently, 
the prediction of thermodynamic quantities for a fluid in this region turns out to be wrong. The accuracy 
of SPT2b also decreases when fluid and matrix particles are of comparable size jl7ll . 

In the current report we improve the SPT2 approach and propose new approximations extending 
them to the case of a one-dimensional system of a hard rod (HR) fluid in a hard rod matrix. It should be 
mentioned that in a bulk one-dimensional HR fluid, the SPT theory |2l|] reproduces the exact result by 
Tonks 1 22] . Some analytical results for thermodynamical properties of a HR fluid in a disordered hard rod 
matrix obtained by the SPT2 approach have already been presented by us in 12311 . However, these results 
have not been tested by a comparison with computer simulation data. That is why to assess different 
approximations in the SPT2 theory we performed grand-canonical Monte Carlo (GCMC) simulations. We 
will show that for a correct description of a HR fluid in a HR matrix it is important to introduce the 
third porosity <p* defined by the maximum possible value of packing fraction of a fluid confined in a 
matrix. It is observed that the new approximations lead to much better results than SPT2b. The obtained 
results will be compared with the expressions for the chemical potential of a HR fluid in a hard rod (HR) 
and overlapping hard rod (OHR) matrices obtained by Reich and Schmidt using the density functional 
theory Q. 

It is worth noting that contrary to the three-dimensional case, the one-dimensional model cannot be 
used for a description of adsorption processes, since there is no flow in such a system due a disconnect- 
edness of pores. Nevertheless, the considered one-dimensional model can be useful in interpreting some 
processes related to a fluid inclusion inside one-dimensional channels such as the transport of water 
and ions through molecular-sized channels in biological membranes 1 25], the transportation of adsor- 
bate molecules through zeolites 1 26], the charge-carrier migration in polymers (27j], etc. An interesting 
example of the realization of one-dimensional model for a fluid was demonstrated in (ijj by confining 
mesoscopic colloidal particles in nonochannels. 

The paper is organized as follows. A brief review of the SPT2 theory for a HR fluid in HR and OHR 
matrices is presented in section 2. In section 3 we formulate the new approximations within the SPT2 
theory, which essentially improve the description of thermodynamic properties of a confined fluid. Some 
computer simulation details will be presented in section 4. In section 5 the comparison of different ap- 
proximations of SPT2 with computer simulation results is presented and discussed. And finally we draw 
some conclusions in the last section. 



2. SPT2 theory for a HR fluid in disordered HR and OHR matrices 



The basic idea of SPT is to insert an additional scaled particle of a variable size into a fluid. This 
procedure is equivalent to a creation of cavity which is free of any other fluid particles. The central point 
of the SPT theory described in jlj, llil - lilll is a calculation of the excess chemical potential of a scaled 
particle, fif 1 , which is equal to a work needed to create the corresponding cavity. In the presence of porous 
medium the expression for the excess chemical potential of a small scaled particle in a one-dimensional 
HR fluid is as foUows S: 



/^s-ln(piAi) = 
-lnpoCRs)-ln 



Rs 

p {Rs)-r)\l+ — 



-In 

l/ld + fla/Bl) 



Po(Rs) 



(2.1) 
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where /3 = 1/kT, k is the Boltzmann constant, T is a temperature, R s is the "radius" (half of length) of 
scaled particle, R\ is the "radius" of fluid particles, 771 = 2pii?i is the fluid packing fraction, pi is the fluid 
density, Ai is the fluid thermal wavelength. The term po(R s ) = e~^ s is defined by the excess chemical 
potential pS? of the scaled particle confined in an empty matrix (limit of infinite dilution) and it has the 
meaning of the probability to find a cavity created by this particle in the matrix in the absence of fluid 
particles. The term pi/o(-R s ) = 1 - 771 (1 + RsIR^Pq 1 ^) is the probability of finding a cavity created by 
the scaled particle in the fluid-matrix system under the condition that the cavity is located entirely inside 
a pore occupied by the scaled particle. 

For a large scaled particle {R s 3= 0) the excess chemical potential p|r x is given by a thermodynamic 
expression for the work needed to create a macroscopic cavity inside a fluid and it can be presented in 
the following way 

2.R 

p>f=u/(i? s )+/3P— (2.2) 

where P is the pressure of a fluid. It is worth noting that this expression has the same form as in a 
bulk case. The difference is only in the presence of a multiplier in the last term l/po(.R s ). This multiplier 
appears due to an excluded volume occupied by matrix particles. For the bulk case this multiplier is 
absent, since po(.R s ) = 1. However, it was also erroneously skipped in the SPT1 theory 1 14] leading to 
the inconsistency discussed and corrected in jl7ll . As it was already mentioned, po{R s ) is a probability of 
finding a cavity created by the particle of radius R s 3* in the matrix in the absence of fluid particles. This 
probability is directly related to different types of porosities of a matrix. First, in R s = R\ 

p (R s = R 1 ) = e-^= ( p (2.3) 

one obtains the probe particle porosity, which is defined by the excess chemical potential of fluid particles 
in the limit of infinite dilution Thus it depends on the nature of the fluid under study. The second one 
is the geometrical porosity, which depends only on the structure of a matrix and is related to the volume 
of the void between matrix particles, and it is defined as po{R s ) in R s = 0: 

p o (i? s = O)=0 o . (2.4) 

The corresponding expressions for po (R s ) are derived under condition R s =£ for the case of HR matrix 

polRs) = l-??o(l + -||) (2.S) 

and OHR matrix 

po(#s) = exp -wjl + ljj) - (2-6) 

respectively, where 770 = 2poRo, po = No IV, No is the number of matrix particles, V is the "volume" 
(i.e., the length in one dimension) of a system, Ro is the "radius" (half of the length) of a matrix particle. 
Therefore, the geometrical porosity for a HR matrix has the form 

<£o = PoU?s = 0) = I-770 (2.7) 

and for a OHR matrix it is 

0o = Po(#s = O)=e- ,?o . (2.8) 
For the probe particle porosity, cf>, the following expression is derived 



(p= (l-77o)exp 
in the case of HR matrix, and 



I-170 



(2.9) 



= exp [-770(1 +t)] (2.10) 

in the case of OHR matrix, where t = RiIRq. As one can see (p is always less than (po, except the limit of 
t = 0, in which <p = <p$. 
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According to the ansatz of SPT, w s [R s ) can be found from the continuity of both expressions <2.1> and 
T2S for pf in R s = 0: 

m 



w{R s ) : 



-In 1 



0o 



(2.11) 



After setting R s = the expression i2.2) yields the relation between the pressure, P, and the excess 
chemical potential, p^ x , of a fluid in a matrix: 



Pi 



Using the Gibbs-Duhem equation 

I dP 
\dp\) T 

an expression for compressibility is obtained 



Pi 



dpi 
dpi i T 



P 



dP 



1 



771 /0o 



dpi) T l-77i/0 (l-77i/0)(l-77i/0 o ) ' 



(2.12) 



(2.13) 



(2.14) 



which makes it possible to obtain the total chemical potential, /3pi = In(piAi) + /3/ij X , after the corre- 
sponding integration of J2.14> over the fluid density. Therefore, as it was shown previously in [23], the 
expression for the excess chemical potential as well as for the pressure can be derived 



. l-J7i/0 



p>f- Ml ) = -lnll-^ 



j 0o-0 l-J7i/0o 



Pi 



00 , l-77l/0 



(2.15) 



(2.16) 



771 00-0 l-77l/0o 

From H2.151 the virial (density) expansion for the chemical potential can be easily obtained in the follow- 
ing form 
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From this expansion it is seen that the second virial coefficient is 

( 1 1 
B z = 2R 1 - + — 

10 00 



(2.17) 



(2.18) 



Later it will be shown that the second virial coefficient obtained from the SPT2 approach is in a good 
agreement with computer simulation results. However, the higher virial coefficients are strongly overes- 
timated. This problem is caused by an assumption that the multiplier l/po (i? s ) introduced by us in the 
expression d2.2t does not depend on the presence of a fluid. Due to this in the considered approach, the 
correct value is obtained only for the second virial coefficient. 



3. Some corrections and improvements in SPT2 theory 

In order to improve the SPT2 approach we consider two approximations proposed in 1 17]. In the first 
of them called SPT2a the porosity in <2.12> is replaced by 0o. As a result one obtains 

l 00/ l-7?l/0O 

pSPT2a ^ 

p = (3.2) 

Pl l-77i/0o 
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The SPT2a approximation leads to the following density expansion for the chemical potential of a HR 
fluid 

00 2 0g 3^ 



(3.3) 



from which it is seen that SPT2a does not reproduce the second virial coefficient given in d2.18t . 

The second approximation called SPT2b can be derived by replacing the logarithm in the second term 
of d2.15t with the linear term of its expansion 







00 ~0 



■In 
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(3.4) 



Therefore, the expressions for the chemical potential and for the pressure of a HR fluid in a matrix can 
be rewritten as follows: 
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(3.5) 
(3.6) 



From the SPT2b approximation the density expansion can be obtained 
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(3.7) 



As is seen, the SPT2b approximation gives the same second virial coefficient as in the common SPT2 ap- 
proach. However, the higher virial coefficients in SPT2b are smaller than in SPT2. Therefore, we can 
expect that, similar to the three-dimensional case of a hard sphere fluid in a matrix 1 17], the SPT2b ap- 
proximation is better than the original SPT2. 

Thermodynamic quantities obtained in the SPT2 and SPT2b approximations have divergences at high 
fluid densities. From expressions <2.15) - t2~16l and <3.5) - t3T6l it is seen that two divergences appear at 
r]i = and 771 = 0o, respectively. Since < 0o, the first divergence (771 = 0) occurs at lower densities 
than the second one. That is why overestimation for the chemical potential in the SPT2b approximation 
appears very quickly when the region of high densities is reached. To avoid this divergence, the first 
logarithm in the expression for chemical potential d2.15> (or d3.5» is expanded as well 



-In 1 
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(3.8) 



Therefore, a new approximation is obtained and hereafter it is called SPT2bl (the first correction for 
SPT2b). The chemical potential and pressure of a HR fluid in a matrix in the SPT2bl approximation are 
as follows: 
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From the expansion of chemical potential in SPT2M 

at ex n>SPT2bl ( 1 1 } ? 



20Q 00o , 



171 
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(3.9) 



(3.10) 



(3.11) 



it is clearly seen that this new approximation gives the same second virial coefficient as in SPT2 and 
SPT2b. The next virial coefficients are smaller than in SPT2 and SPT2b and one can expect that the SPT2bl 
approximation is better than other previous approximations. 
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It should be noted that the expressions for a chemical potential of a HR fluid in HR and OHR matrices 
similar to d3.9> were also obtained by Reich and Schmidt 1 24] using the density functional theory If we 
rewrite the expressions of Reich and Schmidt (RS) in terms of the porosities and 0o, they will have the 
form, which is more general and does not depend on the type of a matrix (HR or OHR) 

11 i; I 0oJ l-77i/0o l-i7i/0o 

One can see that (3.121 differs from (3.9) only by the coefficient in the last term of this expression. As is 
seen in the limit — ► 0o, when 

to *0 =ta [l + *>_l)„*!_l f (3.13) 

I J 

the expression 13.121 exactly reduces to J3.9L but it does not reproduce the second virial coefficient ob- 
tained by the SPT2 approach f2TT8> . It is seen that from (3T9) and (TW\ the new approximation (SPT2M) 
does not have divergence in 771 =0, while the divergence in 771 = 0o remains. 

Also, it is worth noting that from the physical (or rather geometrical) point of view the divergence 
should appear near the maximum value of fluid packing fraction, 77™ ax , available for a fluid in the matrix. 
It is well known that for a bulk one-dimensional system the maximum packing fraction of a fluid is equal 
to one, rff 3 * = 1 (2^)2^]. When a fluid is confined in a matrix this value should be smaller 

rj? iax = 0*, 0<0*<0 O . (3.14) 

For a one-dimensional fluid in a HR or OHR matrix the exact expressions for 0* were obtained by Reich 
and Schmidt (3. In terms of the porosities 0o and 0, the porosity 0* can be presented in a general form 
as follows: 

. = i„ «». (3.15) 

00-0 

This expression can be rewritten in a more compact form 

* [V-Vq] 



v — , 

In vl vq 

where we introduce the quantities, which are inverse to the porosities 



(3.16) 



1 1,1 
y = — , v = — , v = — . (3.17) 

00 0* 

At the same time the quantity v* corresponds to the minimum of specific volume for one fluid particle 
under a close packing condition. In the same way we can interpret vq and v, but for the fluid densities, in 
which the corresponding divergences appear for the chemical potential and pressure, i.e. in r\\ — 0o and 
r\\ = 0. The logarithm in 13.16L which can be presented as 

In — = - In = p [m [R s = Ri) - m {R s = 0)] (3.18) 

Vq 00 

is equal to the difference between the work needed to insert a scaled particle of radius R s = R\ (fluid 
particle) and the work to insert a scaled particle of radius R s = (point particle). To introduce 0* into the 
expression for the chemical potential of a fluid <3.5t the first term is modified in the following way 

?7l(0*-0) 



0*0(l-*7l/0*) 



0*0(l-?7l/0*) 

Such a procedure leads to a new approximation which hereafter is called SPT2b2 (the second correction 
to SPT2b). Therefore, the corresponding expressions for the excess chemical potential and pressure are 



23607-6 



A hard rod fluid in a disordered porous medium 



as follows: 
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(3.21) 



The density expansion of the chemical potential of a fluid in the SPT2b2 approximation can be derived as 



rnSPT2b2 



1 1 

00 



1 1 0*-0 

20* 2 0| 0*20 



1 1 

30* 3 -f" 3 



3 0* 3 



(3.22) 



The third correction for SPT2b is obtained by expanding the first logarithm in <3.20> and it is referred 
to in this paper as the SPT2b3 approximation. This procedure is similar to that in <3.19> . thus the excess 
chemical potential and pressure of a fluid in a matrix in the SPT2b3 approximation can be easily derived 
as follows: 
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(3.24) 



The SPT2b3 approximation yields the corresponding density expansion of the chemical potential 
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(3.25) 



From f3T22> and (3T25) it is seen that the both approximations SPT2b2 and SPT2b3 reproduce the second 
virial coefficient presented in (2.181 . The principal difference of the last two approximations compara- 
tively with the other approximations considered in this study is that the SPT2b2 and SPT2b3 approxi- 
mations have the divergence in the correct value of packing fraction 771 = 0*, which corresponds to the 
maximum value of packing fraction of a HR fluid in a matrix. Therefore, one can expect that these ap- 
proximations provide a correct description of a HR fluid in a matrix at high fluid densities. The accuracy 
of all the considered approximations will be tested by a comparison with computer simulation results. 



4. Computer simulations details 

The grand canonical Monte Carlo simulations (iill were performed within this study in order to verify 
the accuracy of the theory proposed. To this end, a one-dimensional system of hard rods was considered. 
The system consisted of two components. The first component is described by unmovable hard rod par- 
ticles of a matrix. The second one is described by hard rod particles of a fluid that could move in a space 
not occupied by the matrix particles. A hard rod interaction between particles in the system is defined by 
radii Rq and R\, which are half sizes of matrix and fluid particles respectively. We considered two types 
of a matrix with different structures. A structure of the first matrix is modeled as an equilibrated one- 
dimensional HR fluid, thus it is referred to as a HR matrix. The second one presents HR particles which 
are randomly distributed in space and can overlap, and this model is referred to as an overlapping hard 
rod (OHR) matrix. 

To compose configurations of a HR matrix, the preliminary simulations were carried out in canonical 
ensemble. For this purpose, a number of matrix particles equal to Nq = 10000 were taken and after 
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equilibration had been reached the coordinates of them were stored to use them in the main simulations 
for a confined fluid. In the case of a OHR matrix, the coordinates of matrix particles were randomly set. 

The different sizes of matrix particles were considered in simulations, starting with the limit case 
of zero size and finishing with the large enough matrix particle sizes (i?o = O.Ocr; Ro = 0.5a; Rq = 1.5a; 
Ro = 2.5a). It should be noted that all sizes and lengths in our study are presented in units of a, which is 
equal to a diameter of fluid particles. Thus, the radius of fluid particles is i?i = 0.5a. Another important 
parameter for a matrix is the porosity <p, which was fixed in our study and was equal to 0.35. Therefore, 
depending on the size of matrix particles, the density of matrix particles was adjusted to make <p equal to 
this value using the expressions (279) and <2.10> for HR and OHR matrices, respectively. Since, a number 
of matrix particles, No, were fixed, the volume of the box depended on the matrix density: 

N 

V = L=— , (4.1) 
Po 

where po-rjoi 2Rq . A number of fluid particles varied depending on the chemical potential of a fluid up 
to Ni ~ 15000. During each simulation run, equilibration and production procedures were performed. 
The number of steps for equilibration equal to 600000, and 200000 steps for the production appeared 
to be sufficient to get reliable results for the fluid density at the given chemical potentials. It should 
be also noted that the results show a remarkable dependence on a matrix configuration. That is why 
in such a case, different configurations of a matrix should be used to take an average over them. To 
this end, we used 8 matrix configurations for each set of matrix parameters and the chemical potential 
to get a satisfactory statistical error for the fluid density, which was less than 0.5%. Therefore, from the 
simulations we obtained numerical relations between the fluid density and the chemical potential, which 
were compared with the results calculated from the SPT2 approach. The corresponding analysis is made 
in the next section. 



5. Results and discussions 

Above we presented the analytical expressions for the chemical potential and the pressure, which 
were obtained in different approximations for a one-dimensional hard rod (HR) fluid confined in a dis- 
ordered matrix. These expressions include three types of porosities, i.e. the geometrical porosity (po, the 
specific probe particle porosity <p and the porosity <p* defined by the maximum packing fraction of a fluid 
in a given matrix. All these porosities (po, <p and <p* have exact expressions I2.7> — f27T0t and I3.15t . 

Before testing the accuracies of the proposed approximations we consider two limits. One of them 
t = R\IRq — ► was discussed previously in fl7ll for a hard sphere fluid in a matrix. In this limit, a size 
of matrix particles tends to infinity, Ro — ► oo, while the packing fraction of matrix particles 770 is fixed. In 
this limit one can see that (p Q -(p-(p* and the considered approximations lead to the same result, which 
corresponds to a bulk fluid with a density 771 = rj\/(po. In the second limit t — » 00 and a size of matrix 
particles tends to zero, Rq — * (770 — ► 0), while 770 t = 2poR\ is kept constant. This limit corresponds to the 
matrix of point particles and in this case the difference between HR and OHR matrices disappears. The 
expressions of porosities for a point matrix are as follows: 

0o = l, cf, = e -»°\ ^ = -^-r. (S.D 

Now we proceed to assess the accuracy of different approximations of SPT2 at different values of T 
and fluid densities. For convenience, we fix the probe particle porosity, cf>, and for all cases considered in 
this study we set (p = 0.35. For t we consider three cases t = 1, r = 1/3 and t = 1/5. We start with a com- 
parison of the approximations SPT2, SPT2a and SPT2b with the results of grand canonical Monte-Carlo 
simulations obtained in this study (figure 1). Similar to the three-dimensional case 1 17] SPT2a underes- 
timates and SPT2 overestimates the values of chemical potential in comparison with GCMC. The SPT2b 
approximation improves the results for low and intermediate fluid densities. With a decreasing t the 
improvement becomes more essential. Since SPT2 and SPT2b have the divergence at 771 = </>, the descrip- 
tion is not correct at higher densities. This drawback should be eliminated by the new approximations 
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SPT2M, SPT2b2 and SPT2b3. In figure 2 a comparison of these approximations with the SPT2b approx- 
imation and with the expressions by Reich and Schmidt [24] are presented (in the figure denoted as RS 
results). As one can see SPT2M essentially improves the results at intermediate densities, but with a re- 
markable underestimation at high fluid densities. It is observed that the RS results are rather close to that 
obtained in SPT2bl, but it yields a little bit lower values of the chemical potential. With a decrease of t, 
the difference between these two approximations reduces. However, due to the divergence in 771 - (po 
appearing in these two approximations, the chemical potential can be finite at r\\ > (p* {(p* < (po) and 




Figure 1. The excess chemical potential of a one-dimensional HR fluid in a disordered HR (left panel) or 
OHR (right panel) matrix. A comparison of SPT2, SPT2a and SPT2b approximations (lines) with the GCMC 
simulation results (symbols). 
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this is not correct since the packing fraction of a fluid cannot reach the values that are higher than the 
maximum packing fraction available for a fluid in a matrix, (p* . The best description is provided by the 
SPT2b2 and SPT2b3 approximations, which have a correct divergence at rji = <p* , while SPT2b3 is slightly 
better than SPT2b2. 





Figure 2. The excess chemical potential of a one-dimensional HR fluid in a disordered HR (left panel) or 
OHR (right panel) matrix. A comparison of different approximations (lines) with the GCMC simulation 
results (symbols). The vertical dotted line corresponds to the value of maximum packing fraction of a 
one-dimensional fluid confined in a HR or OHR matrix, (p* . 



The results for a one-dimensional HR fluid in a point matrix (t — oo) are presented in figure 3. As one 
can see, in this case, SPT2b is correct only up to the density 771 = 0.2. SPT2bl provides a correct description 
up to 771 = 0.45. In this case, the difference between SPT2bl and RS approximations is more pronounced 
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and SPT2bl shows a better agreement with the simulations, while SPT2b2 and SPT2b3 provides a the 
correct description up to 77 1 = 0.3. They also provide a qualitatively correct description of the chemical 
potential at high fluid densities. 
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Figure 3. The excess chemical potential of a one-dimensional HR fluid in a disordered point matrix. A 
comparison of different approximations (lines) with the GCMC simulation results (symbols). The vertical 
dotted line corresponds to the value of maximum packing fraction of a one-dimensional fluid confined in 
a point matrix, <p* . 




1.00 \ 1 , 1 , 1 , 1 , 

0.00 0.01 0.02 0.03 0.04 

*1, 



Figure 4. The excess chemical potential of a one-dimensional HR fluid in a OHR matrix at very low den- 
sities. Solid line denotes the SPT2 approach and dotted line corresponds to the result obtained by the RS 
formula (24b . 



Finally, we consider the description of chemical potential at small densities, where it is sufficient to 
take into account only the second virial coefficient. In sections 2 and 3 it is shown that for all of the con- 
sidered approximations, except the SPT2a and RS, the second virial coefficient has the same expression 
d2.18t . In figure 4 the dependence of the chemical potential of a HR fluid in a OHR matrix on a fluid den- 
sity is shown at the small densities up to r\ — 0.04. As one can see, in this region the chemical potential 
is correctly described by the second virial coefficient in the form d2.18t . The second virial coefficient ob- 
tained from the RS approach leads to a slight underestimation of the chemical potential in the region of 
the considered fluid densities. 
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6. Conclusions 

In this paper the recently proposed scaled particle theory SPT2 1 17] for the description of hard-sphere 
fluid in disordered porous media is applied to a one-dimensional hard rod fluid in HR and OHR matrices. 
The obtained analytical expressions for the chemical potential and pressure of a HR fluid in a disordered 
matrix explicitly include two types of porosity. One of them is geometrical porosity, (po, describing the free 
void in a porous media, and the second one is the probe particle porosity, <p, described by an adjustable 
volume for a fluid in a porous medium and is defined by the chemical potential of a fluid in the limit of 
infinite dilution. It is shown that the SPT2 approach correctly reproduces the second virial coefficient, 
but overestimates the higher virial coefficients. In order to improve the SPT2 theory we have applied two 
approximations - SPT2a and SPT2b. To assess these approximations, the extensive GCMC simulations 
were carried out. A detailed comparison of the results of the approximations with the simulations shows 
that all of them have got a good accuracy at small and intermediate fluid densities. However, SPT2a leads 
to a remarkable underestimation of the chemical potential while SPT2b causes overestimation at high 
densities. 

The principal defect of the SPT2 approach is related to the divergence of thermodynamical properties 
at the densities corresponding to the probe particle porosity r\\ — (p. In order to avoid this defect three 
new approximations were proposed. In one of them called SPT2bl, thermodynamical properties have a 
divergence only at the density corresponding to the geometrical porosity 771 = <po. This approximation 
provides a very accurate description at rather high densities for different parameters of porous media 
including the matrices formed by disordered point particles. To describe thermodynamical properties in 
the region of very high fluid densities we have formulated two new approximations SPT2b2 and SPT2b3, 
which explicitly include the third porosity <p* defined by the maximum value of packing fraction of a fluid 
in a matrix. These approximations have a divergence at 771 = <p* , which provides a correct description of 
thermodynamical properties at high densities. 

It is worth noting that the theory proposed and validated in this report can be easily extended to the 
cases of two- and three- dimensional systems. The application of the developed theory in combination 
with the new approximations for a three-dimensional hard sphere fluid in hard sphere matrices will be 
described elsewhere. 
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OflHOBMMipHMM TBepflOCTpM>KHeBMM 11/11/1 H y 

HeBnopfiflKOBaHOMy nopi/icTOMy cepeflOBi/imi: 

Teopifl MdCLLITd6HOI MaCTMHKM 

M. To/iobkcP, T. naL^arai-P, B. flom 

iHCTHTyT (J>i3HKH KOHfleHCOBaHWX CWCTeM HAH YKpaiHM, YKpaiHa, 79011 /lbBiB, By/1. CBCHL4iL4bKOrO 1. 

Bi/ima HopMa/ibHa WKO/ia JliOHy, /la6opaTopia xiMii, <t>paHU,ifl 

3acTOcoBaHO Teopiio MacujTa6H0T nacTUHKH flo onucy TepMOflUHaMmHwx B/iacTHBOcreii n/iwHy b MaTpwuj b 
paMKax MOfle/ii TBepfli/ix CTepxmB. 3 u,i£K) MeTOio po3BHHyTO TeopeTH4Hwti niflxifl, BiflOMnPi ak SPT2. OipwMa- 
ho aHa/iiTH4Hi Bwpa3w fl/ia xiMi^Horo noTeHuja/iy i TUCKy n/iwHy b MaTpnuj flBOx TuniB: TBepfli crepxHi Ta TBepfli 
CTepxHi, mo nepeKpi/iBatOTbCH. 3anponoHOBaHO HH3Ky hobwx Ha6/inxeHb b paMKax SPT2. riOKa3aHO, mp, k p i m 
BiflOMWx reoMeTpn4HOi nopucTOCTi Ta nopwcrocTi npo6HOi 4acTHHKH, Baxm/iBWM e me iHwuw ™n nopucrocTi, 
akhw BW3Ha4aeTbCfl MaKCHMa/ibHO mox/iubok) ynaKOBKOK) n/iwHy, flocrynHOK) b nopncTOMy cepeflOBumj. LHo6 
nepeBipnTH Teopiio SPT2 Ta oujHw™ TOHHiCTb 3anponoHOBaHnx Ha6/inxeHb, npoBefleHO KOMn'KvrepHe mo- 
,qe/iK)BaHHfl MeTOflOM MoHTe-Kap/10 y Be/WKOMy KaHOHHHOMy aHcaM6/ii. 3ayBa>KeHO, 140 TeopeTUHHUii onnc, 
npeflCTaB^eHww yflaHOMy flocniflxem-ii, cyrreBO noKpamye pe3y/ibTa™ ax,qo HafiBWLUnx rycTWH rmi/iHy. 

K/iK)HOBi c/ioBa: n/inHH s npocropoBfix o6Me>K.eHHRX, nopucri Marepia/in, reopia MacwTadHO'i nacimm, 

TBepflOCTptDKHeBHW rWI/IH, TepMOflUHaMi^Hi BJiaCTWBOCTi, KOMn'K)TepHe M0flejlK}B3HHfl 
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